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SEM STUDIES OF WHITLOCKITE AND CHROMITE IN METAL- SULFIDE DROPLETS IN 
A C3 ALLAN HILLS METEORITE AND THEIR IMPLICATIONS. Eun-Hee Cirlin. Rockwell 
International Science Center, Thousand Oaks, CA 91360. 

The meteorite ALHA 77003 is classified as type C3(l) and was referred to 
as type C03(2). In this study whitlockite and chromite were found in metallic 
droplets in both chondrules and in matrix of polished section 77003,78 and were 
studied with a scanning electron microscope. Our SEM is equipped with a Prince- 
ton Ganna Tech x-ray analyzer, which is Interfaced with a Kevex 7000 data anal- 
ysis system. The quantitative analyses with standards using the program MAGIC 
V give better than 2% accuracy for major elements and satisfactory results for 
minor elements. 

Recently Rambaldi et al (3) reported finding whitlockite in unequilibriated 
ordinary chondrites for the first time. They found the whitlockite in the 
interfaces between Si, Cr-rich kamacite and sulfides. In the present SEM stud- 
ies, the whitlockite was found in all three phases, kamacite, taenite and troi- 
lite of metal-sulfide droplets in both chondrules and matrix. Many whitlockite 
worm-like grains were present along the edges of the metallic droplets (Fig. 
la). Troilite was present either as a shell around the metal (Fig. la), or with 
sharp crystal boundary lines (Fig. 1c). Also other small Inclusions (< 3 nm) 
of chromite and Fe-rich olivine were present in the metal (> r ig. lb and 1 c) . 

Table 1 shows analysis of typical whitlockite. 

It was suggested for Apollo 14 samples that whitlockite was formed as a 
result of a reaction of schreibersite with surrounding pyroxene and oxygen (4). 
Our SEM studies showed that Ca in chondrule pyroxene, next to the whitlockite 
tooth in troilite is indeed depleted (Fig. Id), suggesting that a similar reac- 
tior occurred in the meteorite. Since phosphorous has little solubility in sul- 
fides (5), it must have been in the metal before the whitlockite grains were 
formed by oxidation. Then the metal was partially sulfurized some time later. 
Therefore, finding of similar whitlockite grains in both metal ar.d trol 1 i te 
imply that the troilite formed from metal in a later stage after the whitlockite 
grains were formed and therefore after the chondrules were formed. Theoretical 
calculations by Olsen and Fuchs (6) supported by experimental work of Friel and 
Goldstein (7) can be used to show that in thv; solar nebular condensation pro- 
cess, the whitlockite would have formed at less than 680 ft K, whereas the troi- 
lite would be formed at around 70C°K. This suggests that the troilite in ALHA 
77003 was not formed during the solar nebular condensation sequence but was pro- 
bably formed at a later time in a planetesimal (8). 

The conclusion from the present study of this unequilibriated meteorite is 
that the *roilite in chondrules resulted from sulfurization which occurred after 
the chondrules formed. This is in strong contrast to the conventional interpre- 
tation (9, page 179). 

Table 1. Energy Dispersive X-ray Analysis of a Whitlockite 
Element Na M^ Si_ £ Ca Fe Total 

No. of cations in formula 07?4 0.32 0.11 1.91 27T6 0.31 5.25 

based on 8 oxygens 

(1) Marvin U. B. and Mason B. (1982) Catalog of Meteorites from Victoria Land , 
Antarctica, 1978-1980, p. 85. (2) Swart P. K. , Grady M. M. , Norris S. J., Wright 

I. P. and PiTlinger C. T. (1983) Lun ar and Pl anet ary Science XIV , p. 763. (3) 

Rambaldi E. R. and Rajan R. S. (198?) Mete oriti cs 1 p. 271. (4) El Goresy A., 
Taylor L. A. and Ramdohr P. (1972) Proc. Lunar Sci . Conf. 3rd. p. 333. (6) Olsen 
E. and Fuchs L. H. (1967) Icarus 6 , p. 2T2. T"5T Keil K. (1968) U. Geophys, Res., 

73, p. 6945. (7) Friel J. J. and Goldstein J. I. (1976) Proc. Lunar Sci. Conf 

7fh, p. 791. (8) Armstrong J. T. , Meeker G. P., Huneke J. C. and Wasserburg G. J. 
TT582) Geochimica et Cosmochimica Acta, 46, p. 575. (9) Dodd R. T. (1981) Meteor- 
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SEM STUDIES OF WHITLOCKITE AND CHROMITE 


ites: A Petrol oqic-Chemi cal Synthesis , Cambridge University press 


Fig. la. Backscattered electron image (BSE) of a metal droplet showing troil- 
ite shell around kamacite. Dark worm-like whitlockite grains are present in 
both kamacite and troilite along the edge of the droplet (b) BSE of small 
inclusions such as chromite and Fe-rich olivine present in kamacite 
(c) SEM photo of a metal droplet with chromite and whitlockite grains. Notea 
sharp boundary between kamacite and troilite (d) A whitlockite tooth in troi- 
lite. Ca in chondrule pyroxene near the whitlockite tooth is very depleted, 
suggesting that Ca from surrounding pyroxene reacted with phosphorous from the 
metal phase forming whitlockite. 
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OXIDIZED, REFRACTORY AND ALKALI-RICH COMPONENTS IN QINGZHF.K 
ENSTATITE CHONDRITE: IMPLICATIONS ABOUT THEIR ORIGIN. E. R. Raabaldl, R. 
S.RaJan, Jet Propulsion Laboratory, California Institute of Technology, 
Pasadena, CA; R. M. Housley, Rockwell International Science Center, Thousand 
Oaks, CA; and D. Wang, Institute of Geochemistry , Academia Slnlca, Gulyang, 
People's Republic of China. 


Recent studies have shown that the highly unequlllbrated enstatlte 
chondrlte Qlngzhen contains a component that appears to have originated under 
relatively high oxidizing nebular conditions (1, 2). This component has been 
only Incompletely reduced and consists of: a) relict grains of dusty olivine 
within Incompletely melted chondrules (1), and b) enstatlte grains, either 
Isolated or in chondrule Interiors (2). In order to Improve our understanding 
of the origin of these materials, we are continuing our detailed chemical and 
textural studies of the silicate portion of Qlngzhen. Some of the results are 
presented below. 

Forsterlte chondrules 


Qlngzhen contains a small population of chondrules which are enstatlte- 
free and consist almost entirely of forsterlte. Chondrule 3 consists of an 
Intergrowth of forsterlte crystals with minor amounts of Ca-pyroxene, 
plagloclase, metal and glass. The glass phase Is small In size (<10 urn) and 
occurs both in interstitial regions and as Inclusions within the olivine. Its 
composition varies from almost pure sodalltlc C 7.5Z Cl), to predominantly 
albltlc, but other compositions having variable Ca and K contents have also 
been found. The Ca-pyroxene forms rims around the various glass phases and 
displays a wide compositional range (A1 2 0 3 5.1-16.2; T10, *2-1.3; CaO 15.5- 
22.9; MgO 13.4-23.5; Cr 2 0. .6-1.6). A few grains of labradorltlc plagloclase 
•re associated with the glass and Ca-pyroxene. This Is the first report of a 
Ca-rlch plagloclase within an E-chondrlte. The metal droplets have low N1 and 
SI and high Cr contents (Ni .03-2.3; SI .12-1.6; Cr .07-6.6). Their 
composition is suggestive of the reduction process of precursor oxidized 
material. The large compositional variation of the glass, metal and Ca- 
pyroxene phases Indicates that. In spite of Its Igneous texture, chondrule 3 
was never completely melted and homogenized. On a SlOj-CaO+AljOj-MgO ternary 
diagram, this chondrule falls Into a region having a higher Mg/sl ratio than 
the enstatlte -rich Qlngzhen chondrules. 

Ca. Al-rlch chondrule 

Chondrule 17 contains euhedral to subhedral enstatlte crystals, 
representing about 1/5 of the chondrule volume, surrounded by a fine-grained 
Intergrowth of labradorltlc plagloclase (An 61-69) and S10 2 -rlch glass. A few 
coarser grains of labradorlte and Ca-pyroxene are scattered throughout the 
Intergrowth. The whole chondrule la enriched In Ca (xlO), T1 (xlO), and A1 
(x5) relative to the enstatlte-rlch chondrules In Qlngzhen. 

In enstatlte chondrites the bulk of the Ca and T1 occurs In sulfide 
phases (oldhamlte and trolllte) and only minor amounts of Ca are present In 
the lnterstltltal glass of chondrulea. This suggests that the precursor 
material of chondrule 17 originated In a leas reducing, leas aulfurlzlng and 
more refractory environment than that which existed at the formation location 
of the enstatlte chondrites. 
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While minor amounts of forsterlte have been reported within some 
chondrulea of EH chondrltea (3), our work on Qlngzhen la the flrat report of 
the presence of forsterltlc chondrulea and chondrulea containing large amounts 
of Ca, Al-slllcates in this type of chondrlte . The presence of these 
materials supports the dust-gas fractionation model suggested by Larimer and 
coworkers (4, 5) as a way to produce reducing conditions at the enstatlte 
chondrlte formation ' location. Thus, the refractory and oxidized components 
found In Qlngzhen might represent materials which originated In dust-enriched 
nebular regions and were subsequently transported into an adjacent, more 
reducing environment, In which the Qlngzhen chondrlte was formed. 


Alblte-roedderlte-metal assemblages 


Qlngzhen contains round metal-sulfide assemblages, some of which appear 
to have experienced reheating and. In some cases, melting. Small amounts of 
SIO 2 and enstatlte often occur as Inclusions within the metal of these 
assemblages. Within sulflde-poor assemblages we have observed substantial 
amounts of alblte and roedderlte Included In the metal phase. In some cases 
the two phases are Intimately lntergrown and small amounts of SIO 2 are also 
present. Neither mineral has been found In other textural locations, with 
the exclusion of albltlc glasses within chondrules. 


As previously observed for Ca .ml Tl, In Qlngzhen the bulk of the alkalis 
occur within major sulfide phases: djerf lsherlte (K), caswellsllverlte (Na) 
and Cr-sulfldes (Na, K) (6, 7). It Is therefore surprising that alkali- 
bearing silicates are found In Qlngzhen In this peculiar textural location. 
It appears that the metallic host has prevented the Na and K-slllcates from 
being sulfurlzed by reaction with nebular H 2 S, thus preserving their original 
compositional and textural characteristics. Even though the origin of these 
alkali-bearing assemblages Is still unclear it Is possible that they may also 
represent components originated In a different, more oxidizing nebular 
environment. 

REFERENCES . I) Rambaldl, E. R. , £t _al. (1984) Earth Planet. Scl. Letters , in 
press. 2) Rambaldl, E. R. , et al. (1984), Nature , submitted. 3) Blnns, R. A. 
(1967) Am. Mineralogist, 52 . p. 1549-1554. 4) Larimer, J. W. and Bartholomay, 
M. (1979) Geochlm. Cosmochim. Acta, 43 . p. 1455-1466. 5) Bartholomay, H. A. 
and Larimer, J. W. (1982) Meteorltlcs. 17. p. 180-181. 6) El Goresy, A. et 
al . (1983), Meteorltlcs . In press. 7) Rambaldl, E. R. ^t^ al. (1983) 
Meteorltlcs . In press. 
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CHARACTERISTICS OF REMELTED ALLENUE CHONDRULES AND MATRIX. 

R. M. Housley and E, H. Clrlln. Rockwell international Science Center* Thous- 
and Oaks, CA 91360. 

As part of a more general observational and experimental Investigation 
Into the conditions under which various components In unequilibrated meteorites 
may have formed we have melted several Allende matrix samples and have remelted 
several Allende chondrules. This was done In the miniature graphite furnace of 
our atomic absorption apparatus (1), which provided conditions as reducing as 
plausible for any chondrule formation environment. Samples were rapidly heated 
to maximum temperatures In the range from 1350* C to 1750* C, held there for .2-3 
seconds* and then cooled at 300-400° C per second. 

The above procedure has the obvious shortcomings that sample compositions 
are not known In advance* and that the cooling rates are not controlled. Des- 
pite these* because of the very realistic nature of the starting materials* It 
provides Important general Information complimentary to that provided by pre- 
vious more comprehensive and well controlled synthetic chondrule experiments 
starting with glass or oxide powders (2-4). It also provides Information di- 
rectly applicable to Interpreting features seen In Allende chondrules. 

The most Important general feature observed In all our runs so far Is that 
Na loss Is almost complete while reduction of Iron Is fairly limited, even 
though the environment Is highly reducing. Although kamaclte Is present In 
most samples, high FeO silicates and pentlandlte always remain. This seems to 
rule out the possibility that chondrules could be produced by melting matrix In 
a reducing environment followed by the separation of reduced metal. 

Except for containing large bubbles of doubtful significance, chondrules 
remelted at about 1500° C generally have textures typical of fresh porphyrltlc 
or ml croporphyrl tl c chondrules as Illustrated In Fig. la. A large fraction of 
the glass In actual Allende chondrules has been altered to secondary phases. 

As expected with this thermal history, most of the forsterltlc olivine In 
these remelted chondrules consists of relatively unchanged relic grains. As 
evidence of this the trapped melt Inclusion shown In Fig. 2b Is rich In Na and 
K, while these elements are undetectable In the nearby mesostasls glass. Little 
grain growth of olivine or crystallization of material from the glass Is evi- 
dent In the remelted chondrule Illustrated In Figs, la and 1b. 

In other cases several microns of olivine grain growth occurred, accompan- 
ied by the growth of micron sized crystals from the melt. This Is Illustrated 
In Figs. 1c and Id. In general we believe the effect of remelting chondrules 
to this degree Is largely to restore them to a conditions closely resembling 
their state before they were subjected to low temperature alteration. 

In these remelted chondrules we never see magnetite, porous pentlandlte, 
or Ca phosphate. These are all things we previously Interpreted to be products 
of low temperature alteration (5). However, we do see virtually unchanged the 
healed cracks containing submicron chromite forcing cs to consider the possi- 
bility that they formed at high temperature, probably by the oxidation of Cr ++ 
which was In the olivine. The presence of 1 urn sized chromite grains as shown 
In 1c and Id, which are never seen In untreated Allende chondrules, Indicates 
that Cr was not in the Cr +++ state when the chondrules formed. 

REFERENCES 

(1) Clrlln E. H. and Housley R. M. (1977) Proc. Lunar Sci. Conf. 8th , p. 3931- 
3940. (2) Blander M. and Katz J. L. (1967) Seochlm. Cosmochlm. Acta 31 . p. 
1025-1034. (3) Klein L. C. , Fasano B. and Hewlns R. ti. (1980) Proc. Lunar 
Plant. Sci. Conf. 11th . p. 865-878. (4) Tsuchlyama A., Nagahara H. and Kushlro 

T (1983) Proc. Confr~Chondrules and their Origins . (5) Housley R. M. and 
Clrlln E. H. (1983) Proc. Conf. Chondrules and their Origins. 
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Fig. 1. Scanning electron micrographs of remelted chondrules. Microporphyri- 
lic texture in one remelted chondrule is illustrated in a and b along with 
submicron chromite along healed cracks. Grain growth of large olivines is 
illustrat d by the backscattered electron Images in c and d. About 1 micron 
of darker, more Mg-rich olivine surrounds the nlivine cores on the right in c 
and on the bottom in d. This in turn is normally zoned to a very Fe-rich rim. 
Bright euhedral 1 urn chromites are evident in the mesostases in c and d. 
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CHARACTERISTICS OF REMELTED ALLENDE CHONDRULES AND MATRIX. 

R. M. Housley and E. H. Clrlln. Rockwell International Science Center, Thous- 
and Oaks, CA 91360. 

As part of a more general observational and experimental Investigation 
Into the conditions under which various components In unequilibrated meteorites 
•nay have formed we have melted several Allende matrix samples and have remelted 
several Allende chondrules. This was done In the miniature graphite furnace of 
our atomic absorption apparatus (1), which provided conditions as reducing as 
plausible for any chondrule formation environment. Samples were rapidly heated 
to maximum temperatures In the range from 1^50* C to 1750* C, held there for 2-3 
seconds* and then cooled at 300-40CP C per second, 

The above procedure has the obvious shortcomings that sample compositions 
are not known In advance, and that the cooling rates are not controlled. Des- 
pite these, because of the very realistic nature of the starting materials, It 
provides Important general Information complimentary to that provided by pre- 
vious more comprehensive and well controlled synthetic chondrule experiments 
starting with glass or oxide powders (2-4). It also provides Information di- 
rectly applicable to Interpreting features seen In Allende chondrules. 

The most Important general feature observed In all our runs so far Is that 
Na loss Is almost complete while reduction of Iron Is fairly limited, even 
though the environment Is highly reducing. Although kamaclte Is present In 
most samples, high FeO silicates and pentlandlte always remain. This seems to 
rule out the possibility that chondrules could be produced by melting matrix In 
a reducing environment followed by the separation of reduced metal. 

Except for containing large bubbles of doubtful significance, chondrules 
remelted at about 1500*C generally have textures typical of fresh porphyrltlc 
or ml croporphyrl tl c chondrules as Illustrated In Fig. la. A large fraction of 
the glass In actual Allende chondrules has been altered to secondary phases. 

As expected with this thermal history, most of the forsteritlc olivine In 
these remelted chondrules consists of relatively unchanged relic grains. As 
evidence of this the trapped melt Inclusion shown In Fig. 2b 1$ rich In Na and 
K, while these elements are undetectable In the nearby mesostasls glass. Little 
grain growth of olivine or crystallization of material from the glass Is evi- 
dent In ths* remelted chondrule Illustrated In Figs, la and 1b. 

In other cases several microns of olivine grain growth occurred, accompan- 
ied by the growth of micron sized crystals from the melt. This Is Illustrated 
In Figs, lc and Id. In general we believe the effect of remelting chondrules 
to this degree Is largely to restore them to a conditions closely resembling 
their state before they were subjected to low temperature alteration. 

In these remelted chondrules we never see magnetite, porous pentlandlte, 
or Ca phosphate. These are all things we previously Interpreted to be products 
of low temperature alteration (5). However, we do see virtually unchanged the 
healed cracks containing submicron chromite forcing es to consider the possi- 
bility that they formed at high temperature, probably by the oxidation of Cr ++ 
which was In the olivine. The presence of 1 urn sized chromite grains as shown 
In 1c and Id, which are never seen In untreated Allende chondrules, indicates 
that Cr was not In the Cr +++ state when the chondrules formed. 

REFERENCES 

(1) Clrlln E. H. and Housley R. M. (1977) Proc. Lunar Scl. Conf. 8th , p. 3931- 
3940. (2) Blander M. and Katz «J. L. (1967) Geochlm. Cosmochlm. Acta 31 . p. 
1025-1034. (3) Klein L. C. , Fasano B. and Hewlns R. H. (i960) Proc. Lunar 
Plant. Scl. Conf. 11th . p. 865-878. (4) Tsuchlyama A., Nagahara H. and Kushlro 

I. (1983) Proc. Confr~Chondrules and their Origins . (5) Housley R. M. and 
Clrlln E. H. (1983) Proc. Conf. chondrules and their Origins. 
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Fig. 1. Scanning electron micrographs cf remelted chondrules. Microporphyri - 
tic texture in one remelted chondrule is illustrated in a and b along with 
submicron chromite along healed cracks. Grain growth of large olivines is 
illustrated by the backscattered electron Images In c and d. About 1 micron 
of darker, more Mg- rich olivine surrounds the olivine cores on the right in c 
and on the bottom in d. This in turn is normally zoned to a very Fe-rich rim. 
Bright euhedral 1 urn chromites are evident in the mesostases in c and d. 
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SCANNING ELECTRON MICROSCOPE STUDIES OF THE VIGARANO CV3 CHONDRITE: 
LARGE FORSTERITE GRAINS. R. M. Housley, Rockwell International Science Center, 
Thousand Oaks, CA 91360. 

Me are pursuing a comprehensive scanning electron microscope SEM study of 
the CV3 chondrlte Vlgarano. Our SEM Is equipped with an energy dispersive 
x-ray detector which Is Interfaced with a Kevex 7000 data analysis system. 

We can obtain satisfactory major and minor element analyses using the program 
MAGIC V. All observations to date are on USNM section 477. 

In general Vlgarano 1$ very complex with most large objects, such as chon- 
drules, showing evidence of one or more episodes of severe shock. On the other 
hand the matrix contains delicate porous structures which could not have sur- 
vived even mild shock. Therefore, as has been shown for Murchison (1), either 
the other components were shocked before being mixed with the present matrix, 
or the matrix Is a secondary alteration product as Is the case for Allende (2). 

During our Initial examination of several porphyrltlc chondrules In the 
section we noticed that some of the larger olivine grains were rounded or 
Irregular In shape whereas most of the smaller ones were euhedral In outline. 
This led us to suspect that only partial melting occurred during the chondrule 
formation process and that these larger grains were relicts, even though they 
were not significantly different In composition from the smaller ones. In 
support of this Interpretation we soon found nearly pure MgAUO^ spinel crys- 
tals near the centers of two of these large olivine grains, although no spinel 
formed In the glassy mesostasls regions. 

Our section contains eight Irregular olivine fragments 600-1700 urn In size 
which have similar composition (Fai >t ,-Fa 3 .o) to the chondrule olivines which 
we believe to be relicts, and which are all larger than the largest euhedral 
chondrule olivine In the section. We are In the process of studying these 
grains In detail In the hope of attaining additional Insight Into the nature 
of chondrule precursor materials. 

Although the Fa number varies some from place to place within these grains 
no well defined zoning trends have been discerned. Calcium ranges from <0.005 
In two grains up to 0.019 atoms per atom of SI. Chromium was detected in six 
of the grains and ranges up to 0.007 atoms per atom SI In two of them. No 
correlations between Fe, Cr and Ca are apparent. Six spinel crystals were 
found In the grain shown In Fig. Id. Four of these can be seen In 1e. 

Although Irregular In outline, five of these grains show a concentration 
of metal and sulfide droplets near their peripheries as Illustrated In Figs, 
la and Id. All of the grains show evidence of alteration and erroslon near 
their borders as Is Illustrated in Fig. If, so the remaining grains may also 
once have had such enrichments. We tentatively Interpret this peripheral metal 
enrichment to be the result of severe shock followed hy annealing. Lobate 
glass Inclusions as seen In Fig. Id may also have been Introduced In the same 
way. This may have been the chondrule forming event. Five of the grains as 
well as many olivines In chondrules also exhibit either parallel fractures or 
mosaic extinction or both as evidence of subsequent milder shocks. 

Careful examination of these grains provides strong evidence that sulfide 
formation occurred after the metal was Incorporated In them. The ratio of 
metal to sulfide always Increases as the border is approached. Droplets In 
unfractured and unaltered olivine never contain sulfides. These must provide 
the best Information about the precursor metal. They Include droplets of N1- 
rlch taenlte, droplets of kamaclte with varying N1 contents and frequently sig- 
nificant Cr and P, and droplets containing both taenlte and kamaelte. The 
olivine grains shown In Fig. la also contains a number of droplets consisting 
of kamaclte plus Nl-rlch schvelberslte as shown In Fig. 1b. In addition. It 
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contains a unique metal grain containing about 60% Cr shown In Fig. lc. Taken 
together these results seem to indicate that metal compositions were highly 
evolved at the time ot incorporation, but that sulfide formation had not yet 
occurred. 

(1) Lambert P. (1933) Lunar and Planetary Science XIV , p. 411-412. 

(2) Housley R.M. and Cirlin E.rl. (1983) Proc. ConT. Chondrules and their Ori- 
gins. 



Fig. 1. Scanning electron micrographs of two large forsteritic olivine crys- 
tal fragments. Peripheral enrichment in metal and sulfides can be seen in 
the whole grain photos a and d. The droplet enclosed by the upper circle in 
a Is kamaclte. The long metal grain enclosed by the lower circle is enlarged 
In c. It contains 60% Cr. Marginal alteration along the lower border of a 
is illustrated In f. Diagonal light streaks are enriched in Cr and Fe. Area 
above long streak Is forsterlte. Dark area Immediately below is enstatite 
which grades shortly In light Fe-rich silicate and then matrix. Three spinel 
grains from the crystal in d are shown near the center of e. 
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STATISTICAL FLAW DETECTION’ 

APPLICATION TO FLAWS BELOW CURVED SURFACES 


R.K.Elsley, K.V.Fertig, J.M.Richardson, R.S.Linebarger 

Rockwell International Science Center 
1049 Casino Dot Riot 
Thousand Oaks, CA 91360 


X. INTRODUCTION 

The detection of the pretence of flaws in structural saterlals 
it the sost laportant function which Non-Destructive Evaluation 
(NDE) performs* At structures are designed to higher performance 
criteria and as safety and life cycle cost factors becoae sore im- 
portant , it be coses necessary to detect smaller and sore difficult 
to find flaws. This paper presents a practical approach to the 
optiaun detection of flaws in the presence of noise signals. A 
decision theoretic approach (described in sore detail in a compan- 
ion paper by Fertlg. ct al. 1 ) is used to derive a detection algo- 
rithm which is adapted to the noise environment in which a parti- 
cular measurement is being made. An automatic procedure for char- 
acterising the noises and developing the optimum detection algo- 
rithm is presented. Two Implementations of this approach have been 
tested on experimental data and show substantial leprovement over 
conventional detection techniques. One Is a flexible algorithm 
used for research purposes, and the other is a real-time algorithm 
suitable for field Implementation. 


This work was sponsored by the Center for Advanced Nondestructive 
Evaluation, operated by the Ames Laboratory, USDOE, for the Defense 
Advanced Research Projects Agency at the Air Force Wright Aeronau- 
tical Laborst^.'les/Materials Laboratories under Contract No. 
W-7405-ENG-S2 with Iowa State University. 
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Ultrasonic measurements of small flaws are usually Halted by 
the presence of one or a nuaber of noise processes. The problea of 
flaw detection (or flaw characterization) is therefore one of dis- 
tinguishing the flaw inforaatlon froa the accoopanylng noise infor- 
aatlon. The statistical approach to NDE Involves taking specific 
account of the statistical nature of the noise processes when de- 
signing flaw detection or flaw characterization algorithas. The 
statistical approach has been used to develop a number of flaw 
characterization algorithms. 2 * 8 This paper describes the use of 
the same approach for flaw detection. 


The first step In developing a statistical approach to flaw 
detection (or characterization) Is to create a measurement model 
which describes both the signals due to the flaw and the signals 
due to noise sources. The seco nd step Is to calibrate the model by 
performing a set of measurements 'diich determine the Instrumental 
properties of the measurement system and the statistical properties 
of the noise mechanisms. The third step is to derive an algorithm 
which performs the detection (or characterization) In an optimal 
manner with respect to the signals and nolseu present. 


The measurement model used in this paper 16 the following 
(expressed in the frequency domain): 


where 


M - XDpAp + XA^j + XAj. + N 


( 1 ) 


M is the measured signal, 

X Is the response of the transducer and Its associated elec- 
tronics, 

Dp Includes the diffraction, attenuation and other factors In- 
volved in the propagation of the sound to the flaw and back, 
Ap Is the scattering amplitude of the flaw, 

A u describes the material noise scatterers, Including propaga- 
tion factors, 

Ap describes "echo*' noise such as geometries or ringdown, and 
N Is random electronic noise. 


This notation differs from that used In Ref. 1 in that Ref. l's 
term "p" is the product of our terms "X" and "D" and that we ex- 
pllciteiy Include a term which describes the "echo noise". 

The derivation of an optlmian detection algorithm Is done in a 
decision theoretic manner. A detailed description of this deriva- 
tion Is given In Ref. 1. A summary follows: The full problem of 

distinguishing all flaws larger than a certain size from smaller 
flaws or no flaw at all generally leads to a mathematically intrac- 
table problem. Therefore, a key simplification has been made which 
leads to detection algorithms which can be Implemented in real time 
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hardware. This assumption is based on the observation that much of 
the energy scattered from flaws comes from localized areas of the 
flat.-, such as the front surface of a volumetric flaw or the tip 
diffraction of a crack. As a result, this scattering Is Impulse- 
like. W have therefore used a simplified flaw model in which the 
Impulse response of the flaw Is modeled as a delta function or a 
set of a small number of delta functions* This leads to detection 
algorithms which are of the form of a convolution of the measured 
signal with a filter function. For the case of modeling the flaw 
as a single delta function, the filter function is: 

<« 

where C(u) is the power spectrum of the noise sources which can't 
be eliminated by prt processing and * indicates complex conjugation. 


Detection then consists of thresholding the quantity: 

S(t) • T M(u) T(w) exp (juit) . (3) 

u 

This corresponds to a type of matched filter for the detection 
of Impulsive scatterers* 


The advantages of the statistcal approach are the following: 

1. More accurate detection and characterization of flaws. 

2. Automatic adaptation to the particular sample being tested and 
to the particular Instrumentation being used. 

3. Statistical algorithms can give confidence measures which can 
indicate the level of confidence which should be placed in any 
given estimate which the algorithm makes. 

Among the practical implications of the statistical approach are: 

1. The ability to detect and size smaller flaws. 

2. The ability to detect and size flaws in noisier materials. 

3. Increased inspection speed because less highly focussed 
transducers are needed to detect a given flaw. 

4. Less exacting requirements on instrumentation, because the 
algorithms adapt to the properties of the instrumentation. 

The limitations of the statistical approach include: 

1. Large flaws can be detected by conventional means. The 
statistical approach is not needed. 

2. Very small flam will not be detected even by the statistical 
approach. There is, therefore, a range of flaw sizes for which 
the technique is appropriate. 

3. The statistical approach requires digital processing of the 
measured data. Currently, this requires more expensive hard- 
ware than analog processing, but the margin is narrowing. 
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III. IMPLEMENTATION OK MEASURED DATA 

In order for ststlstlcal flaw detection tv be of practical use, 
two things wist be proven. First, It wist be shown that a signifi- 
cant Improvement In flaw detection can be achieved on experlment- 
allv measured data, given the limitations of this data. Second, it 
■ust be possible to implement the method in algorithms and hardware 
which operate fast enough to be of use in practical NDE situations. 

In this program, we have developed two implementations of sta- 
tistical flaw detection. The first is a "research algorithm" which 
is written in the ISP signal processing language and provides a 
flexible vehicle for examining the performance of a variety of 
forms of detection algorithms, statistical and conventional. The 
second is a real time algorithm which is implemented in the Digital 
Ultrasonic Instrument (DUI). 


A* Research algorithm 

The first form in which the statistical approach to flaw detec- 
tion has been Implemented is a general purpose algorithm, written 
in ISP, by means of which a number of variations of the algorithm 
were tested and evaluated. The algorithm consists of a setup or 
training phase, followed by a testing phase. 

1. The setup phase consists of the following parts: 

a. Estimate electrical and A/D converter noise from a set of 
waveforms collected at a single location on the sample. The 
algorithm calculates the mean waveform and the variance 
waveform of this set. The variance waveform is used to es- 
timate of the power spectrum of the electrical and A/D con- 
verter noise and can be used to select the amount of signal 
averaging to be used in subsequent measurements so as to 
reduce the amount of these noises to a desired level. 

b. Estimate echo noise and material noise from a set of wave- 
forms collected at a variety of nominally identical flaw 
free locations in the specimen. The algorithm again cal- 
culates the mean and variance waveforms. The mean waveform 
is an estimate of the echo noise (noise which is Independent 
of position). It is saved and subtracted from each new data 

waveform that is acquired. The variance waveform is used to 

estimate of the power spectrisn of the remaining noises end 
is used for C(u) in the detection filter. 

c. The system response X(u) is measured by means of the 

reflection of the sound beam from a flat surface in the far 
field of an unfocussed transducer or at the focus of a 
focussed transducer. 

d. The diffraction D(u) which the sound beam undergoes in 
reaching the flaw location Is calculated using formulae 

developed by Thompson, et.al. 1 * * * * * * * 9 
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Note that if the flaw la sufficiently In the far field of the 
transducer at all frequencies of Interest, then D(u) is the ease 
for flaw and system response measurements and can be ignored. To 
the extent that the flaw Is not in the far field, the diffraction 
can be determined theoretically, as described in the previous para- 
graph, or experimentally, by means of a scatterer of known proper- 
ties at the exact distance of the flaw. Both of these approaches 
are limited because, In practice, It is difficult either to know 
the sound field well enough to reproduce near field effects theo- 
retically, or to have a scatterer whose properties and position are 
known accurately enough to reproduce them theoretically. 

e. The detection filter F(u) (Eq. 2) is now calculated from the 
quantities determined above. 

f. A detection threshold must be set. This can be done either 
theoretically using the measurement model or experimentally 
using reference specimens as is currently done In most NDE 
measurements. Neither approach is entirely satisfactory 
because In neither case can one be sure that the conditions 
In the test piece have been accurately duplicated. Because 
the purpose of our measurements was to study the algorithms 
rather than to detect specific flaws, a threshold was not 
selected, but rather the detection function Itself was 
displayed for analysis. 

2. The test phase consists of the following parts: 

a. A candidate waveform is acquired. 

b. The echo noise Is subtracted from it. 

c. Two detection algorithms are then applied to the data: video 
detection and statistical detection. Video detection is 
performed by rectification followed by low pass filtering 
using a frequency domain Hanning window centered at zero 
frequency and with halfwidth equal to twice the transducer 
center frequency. The statistical detection is performed by 
filtering the signal using the filter F(u) defined above. 

d. A display Is provided of the waveforms at various steps In 
the computation process and of the waveforms output by each 
algorithm. 

B. Real-time algorithm 

In order for a detection technique to be of practical use, it 
must be able to be Implemented in a form which will operate at the 

speed at which conventional ultrasonic inspections are performed. 

In order to demonstrate this capability and in order to provide for 
more rapid testing of the statistical approach, a real-time version 

of these algorithms has been implemented. This part of the work 
was funded by the CANDIS program. 2 * * * * * * * 10 The algorithms were implemen- 
ted on the Digital Ultrasonic Instrument (DUI), 11 which is a high 
speed all-digital instrument for performing sophisticated calcula- 
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tlons an ultrasonic signals. The DU I controls the notion of the 
Ultrasonic Test Bed* 2 in order to scan over the specimens under 
test. The detection algorithns Implemented on the DUI are simpler 
than the research algorithm described above and therefore serve as 
a test of what simplifications can be made in the original algor- 
ithm is order to speed the computations without significantly re- 
ducing the quality of the results. The real-time algorithm also 
consists of a set-up or training phase, and by a testing phase. 

1. In the set-up phase, the DUI guides the operator through the 

series of measurements required to design the detection filter: 

a. The DUI first asks the operator to provide an echo from a 
flat surface in order to determine the system response. The 
DUI provides ar. oscilloscope display of the waveform to 
guide him in his adjustments of the instrument. When the 
operator signals that he has the correct signal, the DUI 
acquires it, calculates its frequency spectrum and inserts 
it into the detection filter. No attempt is made to correct 
for diffraction in this algorithm, as it is assumed that the 
flaw is sufficiently in the far field of an unfocussed 
transducer or sufficiently close to the focus of a focussed 
transducer that such a correction is not needed. 

b. The DUI then asks the operator for a region of the sample 
over which to make a coarsely spaced 6can in order to deter- 
mine the noise present in the sample. It is desirable that 
there be few or no flaws in this region. The DUI scans the 
region, calculates the variance of the noise waveforms and 
thus determines C(u). In this algorithm the echo noise was 
not estimated and saved for the purpose of later subtrac- 
tion, although this capability has been demonstrated in the 
DUI before. 

c. The DUI then calculates the detection filter. 

d. Because the algorithm also calculates the video detection 
waveform, it is necessary to calculate the required lowpass 
filter. In order to do this, the DUI asks the operator the 
center frequency of the transducer, although it could as 
well have measured this from the system response spectrum. 

2. In the testing phase, the DUI performs the following steps: 

a. The Testbed is scanned to the next position to be measured. 

b. The waveform is acquired, with optional signal averaging. 

c. Three detection algorithms are applied to the waveform: 

Peak of the sampled rf waveform 
Peak of the video waveform 
Peak of the statistical waveform 
The three peak values are stored for each point Inspected. 
The sampled rf waveform is Included in the set because it 
would be the easiest to Implement in a digital system. 

d. At the end of the scan, three graphs are plotted. f3ch is 
the peak output of one of the detection algorithms. 
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This section presents two sets of results. First are the re- 
sults of applying the research version of the statistical detection 
algorithm to a set of data in which the level of electrical and A/D 
convert*.. noise have been varied by changing the amount of signal 
averaging used in acquiring the data. Second, we present the re- 
sults of using the real time algorithm while scanning a specimen. 

The specimen used in these measurements is a block of plastic* 3 
in which crack-like flaws have been induced by laser damage (Fig- 
ure 1). The sample has a cylindrically cut surface which simulates 
the borehole of turbine engine components. The flaw on which these 
measurements is based is a crack 0.5 mm in diameter. Optical in- 
spection shows it to be very flat and very nearly circular. The 
normal to the crack is inclined at an angle of 60* with respect to 
the normal to the cylindrical surface. 


A. Results of the research version of the algorithm 

The flaw is fairly easy to detect using a properly selected 
focussed transducer. In order to provide a greater challenge for 
the detection algorithm, the measurements were performed in a le6s 
than optimum manner. First, an unfocussed transducer was used, 
giving a smaller S/N than a focussed one would have. Recall that 
one of the benefits which is expected from the use of the statisti- 
cal approach is that faster scanning will be possible by the use of 
less highly focussed transducers. Second, the transducer diameter 
was large (0.5 in), further decreasing the flaw signal relative to 
the noise. Third, the flaw was located at the first near-field 
null of the transducer at approximately the transducer's center 
frequency (5 MHz). The result is that the flaw signal is only 
approximately 1 LSB in amplitude at the input to the A/D converter. 

Figure 2 shows a set of 6 waveforms acquired with varying 
amounts of signal averaging while the transducer was aimed at the 
flaw. Each is labled with a relative slgnal-to-nolse ratio based 
on the amount of signal averaging used. The flaw is visible at 
34 (is In the higher S/N cases. At the beginning of the waveform is 
some "echo” noise due to the rlngdown of the front surface echo. 

In the lower S/N cases, the quantization noise of the A/D converter 
is clearly visible. 


The analysis of the noise processes gave the following results: 
Figure 3 shows the mean and the power spectrum of the variance of a 
set of waveforms collected with the transducer at a fixed position. 
The meaa contains echo noise, material noise and the flaw. The 
noise power spectrum is a measure of the electrical and A/D conver- 
ter noise. The sharp structure at 12.5 Mhz is due to clock noise 
in the l/D converter (Biomat ion 8100). Figure 4 shows the mean and 
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Fig. 1. Plastic test 
specioen with laser 
Induced crack. 
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Fig* 2* Set of flaw signals with 
varying amounts of noise pro* 
vided by varying the amount of 
signal averaging* 
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Fig. 3. Mean and noise power 
spectrum of a set of waveforms 
collected at a fixed position. 
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Fig. 4. Mean and noise power 
spectrum of a set of waveforms 
collected at various flaw free 
positions. 


the power spectrum of the variance of a set of waveforms collected 
as the transducer was moved to a variety of flaw free locations. 

The mean in this case is the echo noise which can be removed from 
measured data by subtraction. The variance contains all of the 
random noises present in the signal, including material, electrical 
and A/D converter. It is this waveform that is used a6 the basis 
for the noise power epectrim C(u>) in the detection filter. 
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Figure 5 shows the steps involved in determining the system re- 
sponse function. The first two curves ere the reference wsveform 
obtsined from e fist surfece echo end its spectrum. The third 
curve is the celculated diffraction D for the flew location. Note 
the near field null at 5 MHz. The fourth curve is the product XD. 

Figure 6 contains the remainder of the steps Involved in pre 
paring the detection filter. The first curve is the measured noise 
spectrum (from figure 4). Because it is an estimate based on a 
small number of samples, it not as smooth as the expectation of the 
noise process itself. We have therefore smoothed it by means of a 
low pass filter applied to the spectrum (second curve). Finally, 
the filter function F is calculated (third curve). Note that the 
A/D clock noise will be completely eliminated from the measured 
data because the transducer has no energy at 12.5 MHz and therefore 
the filter Ignores this frequency. 

Figure 7 shows the results of applying the algorithm to the 
waveform labeled ”6dB”. The upper curve is the measured waveform. 
The middle curve is the result of applying video (envelope) detec- 
tion and the lower curve is the result of the statistical detection 
algorithm. Video detection does not clearly distinguish the flaw, 
but a significant response due to the front surface ringdown is 
seen. The ststistlcal approach, on the other hand, shows the flaw 
clearly standing out above all the other noises. 
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Fig. 5. Signals used in calibra- 
tion of detection algorithm: 
a) system response from flat sur- 
face echo, b) spectrum |x| of a), 
c) calculated diffraction |d| at 
flaw location, d) |xd|. 


Fig. 6. Additional calibration 
signals: a) Estimated noise 
spectrin C(w), b) smoothed 
noise spectrum, c) filter 
function |(XD)*/C|. 
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Fig. 7. Comparison of envelope 
(video) detection end sta- 
tlstlcal detection for 
S/N - 6 dB. 


Fig. 8. Performance of 
atatletlcal and video 
detection algorlthae 
vs noise level. 

Solid curves are theo- 
retical fits to data. 
Statistical algorithm 
has 16 dB Inherent 
advantage . 
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The results of applying the two algorithms to the six data 
waveforms are summarised In Fig. 8. The horizontal axis Is the 
nominal signal-to-noise of the waveform. The vertical axis Is the 
measured performance of the algorithm, expressed as the peak height 
of the detected signal at the time when the flaw is known to be 
present, divided by the peak height of the detected signal at any 
other time during the measurement. (In computing this measure, the 
detected signal due to front surface rlngdown has bean Ignored, 
since this is o-t present In many measurement situations.) The 
circles show the results for the atatletlcal detection algorithm. 
For the two lowest S/N cases, the flaw was not detected. The 
squares show the results for the video detection algorithm. It Is 
likely that this algorithm did not detect the flaw In any of the 
four lowest S/N cases. The solid curves are a theoretical calcu- 
lation of the expected (S4N)/N for various values of S/N. Each 
curve was fit to the corresponding act of results by varying S 
(l.e. sliding the curve horizontally). The curves for the two 
algorithms are separated by 16 dB, Indicating that the Inherent 
ability of the statistical algorithm to separate signals from noise 
is about 16 dB better than for video detection In this case. 
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Fig. 9 Results of angular scan over 
crack using real tine algorithm. 
Statistical algorithm successfully 
detects crack. 


B. Results of the Real-Time Algorithm 

The real time algorithm was used in conjunction with the Ultra- 
sonic Testbed to scan a transducer over the specimen containing the 
flaw described above. The pivot point of the transducer manipula- 
tor was positioned at the center of curvature of the cylindrical 
surface of the sample and the transducer was then scanned over a 
25* angular range which Included the flaw. Figure 9 shows the peak 
amplitude of the signals produced by each of the three detection 
algorithms as a function of position. Neither the raw digitized 
r.f. vevtform nor the video waveform detected the presence of the 
crack. The statistical algorithm, on the other hand, did detect 
it, with a signal- to-noise of about 12 dB. 

V. CONCLUSIONS 

The method presented here makes use of an explicit knowledge of 
the noise processes in order to design a flaw detection algorithm 
which optimally detects flaws in the presence of such noise. A key 
assumption which makes the approach implementable in a simple form 
is to model the flaw as having an Impulse response function which 
consists of a set of delta functions. 

Experimental results using this statistical approach show a 
significant Improvement In the detectability of a crack- like flaw 
relative to the results obtained with conventional video detection. 
This approach promises to provide a number of advantages In practi- 
cal testing situations, Including detection of smaller flaws, fas- 
ter scanning due to the use of less highly focussed transducers, 
and less need for operator optimization of the measurement process. 


I V I I I 1 1 T I I 

max or nr mavipomm 


j i i I — i > I I J — u.jl ,q, 



i i i r "i r t i i i i i 

na or STATISTICAL 

omcTio* waVmoam 


OS mm 
CAACK 



POSITION IN SOM (OCOAIIS) 


original 

OF POOR 




Rockwell International 

Science Center 

SC507.66MR 


21 


ORIGINAL PAGE 19 
OF POOR QUALITY 


Rockwell International 
Sctone* Center 

SC507.66MR 



REFERENCES •. 

1. K.W. Fertig, J.M. Richardson, R.K. Elsley, Flaw Detection: 
Theory", This proceedings. 

2. K. A. Marsh, J. M. Richardson, R. C. Addison, R. K. Elsley, "The 
Automation of the Born Inversion for Ultrasonic Flaw Sizing", in 
Rev. of Prog, in Quant. Nondestruc. Eval. 2, D. 0. Thompson, D. 

E. Chimenti, Eds. (Plenum, New York, 1983), pp, 975-987. 

3. J.M. Richardson, R.K.Elsley, "Extraction of Low Frequency 
Properties From Scattering Measurements", 1979 Ultra. Sym. Proc. 
IEEE, 79 CHI 4 82-9 , pp. 338-341, (1979) 

4. Dick Elsley, John Richardson, Bob Addison, "Optimum Measurement 
of Broadband Ultrasonic Data”, 1980 Ultra. Sym. Proc. IEEE, 

80CH 1602-2, pp. 916-921, (1980) 

5. J. K. Richardson, R. K. Elsley, "Semi-adaptive Approach to the 
Extraction of Low-frequency Properties from Scattering Measure- 
ments", 1980 Ultra. Sym. Proc., IEEE, pp. 847-851, (1980) 

6. L.A.Ahlberg, R.K.Elsley, L.J. Graham, J.M. Richardson, "Long 
Wavelength Ultrasonic Characterization of Inc n jions in Silicon 
Nitride", Proc. of the DARPA/AFML Rev. of Prog, in Quant. NDE, 
AFWAL-TR-80-4078, pp.656-662, (1980). 

7. R.K.Elsley, L.A.Ahlberg, J.M. Richardson, "Low Frequency Charac- 
terization of Flaws in Ceramics”, Proc. of the DARPA/AFML Rev. 
of Prog, in Quant. NDE, AFWAL-TR-8 1-4080, pp. 151-164, (1981). 

8. R.K.Elsley, "Inversion Algorithms for Crack-like Flaws", in Rev. 
of Prog, in Quant. Nondestruc. Eval., D.O. Thompson, D.E. 
Chimenti, eds.. Plenum, New York, pp. 537-545, (1982) 

9. R. B. Thompson, T. A. Gray, "Analytic Diffraction Corrections to 
Ultrasonic Scatterubg Measurements”, in Rev. of Prog, in Quant. 
Nondestruc. Eval., 2, D. 0. Thompson, D. E. Chimenti, Eds. 
(Plenum, New York, 1983), pp. 567-586. 

10. CANDIS is an acronym for Computer Aided Non- Destructive In- 
spection System. CAND1S is a Technology Modernization Program 
within the Bl-B Contract F33657-81-C-0210. 

11. R.K. Elsley, "The Digital Ultrasonic Instrument", in Rev. of 
Prog, in Quant. Nondestruc. Eval., 2, D.O. Thompson, D.E. 
Chimenti, Eds. (Plenum, New York, 1983), pp. 1487-1500. 

12. The Ultrasonic Testbed is an automated immersion testing system 
developed under Air Force Contract F33615-78-C-5164. 

13. This sample was provided by Dave Hsu of Colorado State Unlv. 

DISCUSSION: 

From the Floort You mentioned getting rid of the echo noise by sub- 
traction. If there's any motion in the transducer, that turns 
into a differentiation that can be a very noisy process. Do 
you anticipate any problems with that? 

R.K. Elsleyi The question was about the practical aspect of imple- 
menting subtraction. Of course, you can do the subtraction 
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cither in the time or frequency domain because it it a linear 
process. 

Now about jiggles in the transducer which produce the effect 
that the signal doesn't arrive at quite the same time when you 
do subtraction. That certainly does happen. Cne of the tech- 
niques thit we have used in some of this processing is to speci- 
fically time align the signal before doing the subtraction. 
Again, it requires a little extra computation capability, but 
by doing time aligning, we can use subtraction in a variety 
of areas such as in this kind of processing and phased array 
imaging. We can often get 25 DB of suppression of unwanted 
signals by careful subtraction. 

P.M. Cammell (Sigma Research): I notice you are looking for signals 

which are below the least significant data recorded. Am 7 
correct that this process actually requires the existence of 
noise and requires that the noise have certain statistics? 

R. K. Elsley: Yes, I would agree. That's not my choice for doing 

such measurements, of course. It was just a trick to use to 
get added signal-to-noise for testing the algorithm, but I 
agree with your statement. 

J . A. Simmons (National Bureau of Standards): Can you use this method 

to improve on time resolution or positioning resolution for 
the flaw? 

R. K. Elsley: Yes. For example, one of the items that is needed 

to do the Born inversion is a very accurate measure of the 
location of the center of the flaw. The correct way to do that 
is by using the low frequency data. We have a statistically 
based approach for making an optimal estimate of where the 
center of the flaw is based on low frequency data. That's one 
that is specifically addressed by this technique. 
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